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1,0 SUMMARY 


An analytical study of the effects of wind tunnel turbulence on turbofan 
rotor noise was carried out. A previously established prediction model was 
used to evaluate wind tunnel turbulence/rotor interaction noise. Tlie predict 
tion model was first extended to include the effects of turbulence anisotropy 
by deriving and incorporating a turbulence spectrum model based on axi- 
symmetric turbulence theory. The extended prediction model was verified 
through extensive data/theory 'Comparisons with scale model and full-scale 
fan measurements, with various inlet turbulence conditions. 

It was found that the classical Ribner/Tucker sudden-contraction theory 
does not adequately describe the turbulence free-stream to fan-face contrac- 
tion effects for the wind tunnel environfflent, and, hence, some uncertainty in 
the prediction of wind tunnel turbulence/rotor interaction noise was identi- 
fied. It was estimated that the predicted levels could be at most 6 dfi too 
low for the highest contractions normally encountered. 

A literature survey of fan inlet turbulence properties was carried out, 
and data ranges were identified for anechoic chambers, outdoor test stands, 
wind tunnels, and atmospheric flight environment. Considerable scatter and 
variability in turbulence properties was found to exist, and "expected 
average" values for the various test site conditions were determined. 

A parametric study of the effects of fan rotor size, blade number, oper- 
ating line, and tip speed on rotor/turbulence noise was carried out. The pre- 
dicted noise levels for the NASA-Ames 40x80 wind tunnel environment were 
compared with both outdoor test stand levels and altitude flyover levels. 

The extended rotor/turbulence noise prediction model was used to predict 
the noise level trends with tip speed, rotor geometry, and environmental 
conditions. 

Results of the study showed the Ames wind tunnel rotor/turbulence noise 
levels to be 15-20 dB lower than outdoor test stand levels for all combinations 
of tip speed, rotor diameter, blade number, and operating line investigated. 

The wind tunnel levels were also compared with flight levels and found to be 


10-15 dB higher, but still lower than the other sources of fan noise (broad- 
band, multiple-pure-tones, rotor/stator, etc.). 

Due to the random variability of atmospheric turbulence properties, it 
was found that a wide range in rotor/turbulence noise levels can be expected 
from different conditions and that the wind-tunnel and flight bands can 
overlap. 

It was concluded that wind tunnel turbulence/rotor interaction noise is 
sufficiently low such that forward flight effects are adequately simulated as 
far as inlet "clean-up" effects are concerned. Although the wind tunnel tur- 
bulence characteristics (small-scale, low intensity) are far different from 
atmospheric flight turbulence characteristics (very large scale, moderate-to- 
high intensity), the net effect on noise is similar. 

It was found that rotor/turbulence noise, for contracting large-scale 
inlet turbulence, is characterized by strong, narrow peaks at blade passing 
frequency and its harmonics, with no significant contribution to the fan 
broadband noise levels. For axial length scales sufficiently large (greater 
than one diameter), the transverse turbulence intensity and transverse length 
scale/blade spacing ratio are the primary parameters which control the rotor/ 
turbulence noise levels for a given fan. 


2.0 INTRODUCTION 


It is now well known that acoustic tasting of turbofan engines in an 
outdoor test stand facility does not yield the same noise characteristics 
as those measured in flight. As discussed in Reference atmospheric 
turbulence is drawn into the inlet during ground static t.isting, undergoing 
substantial elongation and contraction. Because of the large contraction 
of the inflow streamlines, the turbulent eddies convecting with the flow 
appear as long "sausages" as they pass through the rotor. Many blades will 
successively "chop" the same eddy, producing blade-passing frequency tones 
as well as broadband noise. This rotor/turbulence interaction is usually 
strong enough to dominate other sources of fan inlet noise during static 
testing. 

In flight, the fan inlet floO’ c'oes not undergo very much contraction 
because of the aircraft forward motion. The fan inlet contraction ratio, C, 
defined as fan face through-flow velocity divided by flight speed, is not 
too different from unity at normal flight approach speeds. At typical 
approach altitudes (500-1000 ft or 150-300 meters), the atmospheric turbu- 
lence scales or eddy sizes are much larger than on the ground, so that the 
eddies do not appear to the rotor as "sausages" which are successively 
"chopped" by the rotor blades, but may have cross sections as large as, or 
larger than, the inlet itself. Very little rotor/turbulence interaction 
noise is therefore produced in flight, allowing other sources to dominate 
the observed noise spectrum. 

It is desirable to have a "static" or ground-based test facility for 
acoustic evaluation of turbofan engines; one which correctly simulates the 
flight inlet turbulence conditions. Flight testing is not only costly and 
time-consuming, but has severe limits on test conditions, power settings, 
etc., and :.s subject to variability and uncertainty with respect to atmos- 
pheric propagation, ground reflections, and aircraft location. For research 
and development purposes, extensive hardware changes and elaborate instrumen- 
tation and data acquisition equipment are often, required which cannot be easily 
implemented in a flight test. Acoustic testing of full-scale turbofan engines 


in a large wind tunnel ia a viable alternative to flight teating. Flight 
contraction ratioa and aircraft apeeda ^an be aimulated in a wind tunnel* 
under more controlled conditions* with ieat meaaurement capability and flex- 
ibility rivaling that of a conventional outdoor teat stand. 

The NASA-Ames 40x80 wind tunnel facility at Moffett Field, California, 
has recently been used for acoustic evaluation of turbofan engines under 
simulated flight conditions. This wind tunnel facility has turbulence proper- 
ties which are different from both the >U(.uoor test stand and flight altitude 
environment* and it is important to know whether the rotor/turbulence interac- 
tion noise in the tunnel is low enough to provide a true simulation of flight 
noise characteristics. 

The primary objective of the present study was to evaluate the degree to 
which the NASA-Ames 40x80 wind tunnel simulates flight conditions in terms 
of producing sufficiently low levels of rotor/ turbulence interaction noise. 
Guidelines for the ranges Of fan size, blade number, and operating speeds over 
which sufficiently low rotor/turbulence noise levels can be expected were to 
be established. A previously developed theoretical prediction model was to be 
used for performing the analytical study. 

The theoretical model, described in References 2 through 4, utilizes an 
isotropic model of the inlet turbulence, along with sudden-contraction theory. 
Reference 5, to account for inlet flow contraction effects. Experimental 
measurements, however, indicate that the initial turbulence upstream of the 
inlet is not isotropic, e.g.. Reference 6. Additionally, it is desirable to 
be able to specify independently the axial and transverse turbulence intensity 
and length scales. Hence, the theoretical model was to be extended to include 
an axisymmetric initial turbulence spectrum, prior to carrying out the wind 
tunnel evaluation study. 

The scope of the present study consisted of four primary tasks in support 
of achieving the above-stated objective: 

Task 1 - Modification of previously-developed General Electric rotor/ 

turbulence interaction prediction model to include anisotropic 
turbulence. 


4 


Task 2 - 


EsCablithnent: of repreaentativo turbulence properties for three 
acoustic test facilities or conditions: (1) Outdoor Test Stand 

static facility, (2) NASA-Anies 40x80 ft Wind Tunnel, and 
(3) Flight Test. Only existing available data is to be uti- 
lized. 

Task 3 - Evaluation of rotor/ turbulence noise source characteristics of 
a typical high speed fan for each of the facilities/conditions 
listed in Task 2. 

Task 4 - Evaluation of the degree to which the NASA-Ames 40x80 ft 

Wind Tunnel simulates flight conditions from the standpoint of 
acceptably low rotor/ turbulence interaction noise. 

The following sections describe in detail the work performed in the above 
four tasks. A brief review of the rotor/ turbulence noise theory is first 
given, followed by a description of the anisotropic turbulence theory modifi- 
cations. The survey of inlet turbulence properties is then summarized, fol- 
lowed by the parametric study and wind tunnel evaluation. 


ORIGINAL PAGE IS 
OP POOR QUALITY 


‘3.0 REVIEW OF ROTOR/ TURBULENCE NOISE THEORY 


As discussed in References 2 through there are several mechanisms by 
which inlet turbulence might produce noise. The turbulence may be regarded 
as a pattern of vorticity convecting with the flow whose statistical proper- 
ties are known. The nonuniformity of the velocity associated with the tur- 
bulence convecting with the flow produce.? unsteady fluctuations in angle-of- 
atta,i: on the rotor blades, leading to unsteady blade forces and, hence, 
noise radiation. This mechanism is usually referred to as a dipole source 
and is analyzed in detail in Reference 2. 

Further, when the rotor is loaded, i.e., has some steady lift distribu- 
tion, there is a rotor-locked asymmetric flow pattern spinning in the fan 
duct having a fundamental period equal to the blade spacing. Th?ls asymmetric 
pattern itself is an ineffective noise source for subsonic tip speeds, but 
its inteir^<'tion with inflow turbulence leads to a quadrupole source. Insofar 
as Mach number dependence alone goes, the ratio of quadrupole to dipole noise 
should be as M^. The quadrupole source depends on blade loading, whereas the 
dipole source, to first order, does not. The quadrupole mechanism is treated 
in detail in Reference 3. 

Finally, a third mechanism, postulated in Reference 4, arises tihen a 
zero perpendicular velocity boundary condition is imposed at the blade surface 
for the quadrupole model problem. This mechanism is analogous to the noise 
increase experienced when a vibrating tuning fork is brought near a solid sur- 
face. A detailed analysis of this mechanism can be found in Reference 4. 

Based on the above physical picture of noise generation due to turbulence 
incident on a blade row, theoretical analyses were developed (References 2 
through 4) relating the characteristics of the turbulence, design parameters of 
the blade row, and the spectrum of the radiated noise. A spectral representa- 
tion of turbulence which treats the turbulence as a superposition of shear 
waves was used. The turbulence was assumed to be homogeneous and isotropic. 

The effects of inlet contraction were subsequently accounted for (Reference 4) 
by utilizing the sudden-contraction theory of Reference 5, as illustrated in 
Figure 1. 
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A two-dimensional e^v^cade representation of the blade row was employed, 
and the blades were idealized as flat plates. The flow was assumed to enter 
the inlet axially with axial Mach number M 3 , while the cascade translated 
with transverse Mach number M^;, as illustrated in Figure 2. The turbulence 
was treated as a superposition of shear waves of varying wavenumber in the 
axial (x), transverse (y), and spanwise (z) directions. 

The statistical description of the turbulence assumed determines the 
amplitude distribution of the velocity shear waves over the wavenumber spec- 
trum* The X, y, and z components of the shear wave velocities are resolved 
into components normal to the blade chord to determine the unsteady angle- 
of-attack or upwash incident on the rotor blades. 

For the unsteady blade force evaluation, the co^spressible, two- 
dimensional, isolated airfoil lift response functions given in References 7 
and 8 are used. The low-frequency theory of Reference 7 is used for low 
reduced frequencies, while the asymptotic high-frequency theory of Reference 8 
is used for high reduced frequencies*. A correction for airfoil aspect ratio 
is also included, as was suggested in Reference 9. 

As described in Reference 3, the steady blade-to-blade flow field locked 
to the rotor is modeled by representing the blades as a row of translating, 
equally-spaced concentrated point forces. The velocity field due to the 
steady blade loading is coupled with the turbulence velocity field to evaluate 
the resulting Reynolds stresses which are the quadrupole source strengths, 
analogous to the fluctuating Reynolds stresses which form the source strengths 
for jet mixing noise. Since the amplitude of the rotor-locked steady velocity 
field is proportional to the blade loading or Lift coefficient Cj^^, the quadru- 
pole strength is also proportional to and hence so is the associated 
acoustic pressure. 

The General Electric Company computer code for the prediction of rotor/ 
turbulence interaction noise is essentially a modified form of the code pre- 
sented in Reference 3. The incompressible lift response function (Sears* 
function) has been replaced by the compressible response functions described 
above, and the contraction effects as d scussed in Reference 4 have also 
been added. 
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Figure 2. Rotor/Turbulence Interaction Theory Model Geometry and 
Nomenclature. 
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The quadrupole/blade interaction source mechanism (tuning fork amplifi- 
cation effect) has not been included, for several reasons. First, diagnostic 
calculations have shown this mechanism to exhibit "resonances" in its present 
formulation, which do not appear in the experimental data. Second, the ex- 
perimental effect of blade loading (Reference 6) was very small on blade- 
passing frequency, indicating that the dipole source is dominant. Third, the 
formulation in its present form required certain arbitrary modifications in 
order to prevent the predicted noise level from becoming either too large 
or too small, and it was felt that a degree of rigor was sacrificed compared 
to that of the dipole and quadrupole source models. Finally, the computation 
time for the quadrupole/blade interaction contribution was an order-of- 
magnitude larger than that for both dipole and quadrupole contributions com- 
bined. In view of the uncertainties associated with this portion of the for- 
mulation, the additional expense was not considered worthwhile. 

The above-described computer code prediction, prior to the present pro- 
gram modifications, predicts inlet and exhaust radiated power spectra (PWL) 
in proportional bandwidths. Input required is Mg, C , rotor solidity 
(defined as chord c over blade spacing s), inlet turbulence intensity, axial 
length scale, and inlet contraction ratio. The initial turbulence upstream 
of the contraction was assumed to be isotropic. Output includes dipole, 
quadrupole, and total PWL spectra for both inlet and exhaust ducts. 

A sample prediction of rotor/turbulence noise inlet duct PWL spectrum 
and a corresponding measured inlet arc PWL spectrum, taken from Reference 6, 
is shown in Figure 3. The predicted spectrum utilizes the measured axial 
turbulence intensities and axial length scales reported in Reference 6, 
along with a contraction ratio based on the measured ratio of transverse to 

f 

axial intensity and the sudden contraction relationships of Reference 5. It 
can be seen from this comparison that there is much room for improvement in 
the prediction. The basic problem lies in the isotropic turbulence assump- 
tion. The measurements of Reference 6 show a ratio of transverse turbulence 
intensity ut/U^ to axial turbulence intensity Ua/U^ of about 2.0. The 
measured ratio of axial length scale to transverse length scale 
about 300. Using isotropic turbulence followed by a sudden contraction, it is 
predicted that u^/ug^ is approximately the same magnitude as in strong 

contradiction to the measured evidence. 
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1/3-Octave FWL, dB re: 10 watts 








Hence, before the proposed study outlined in the introduction could be 
carried out, the inlet turbulence model had to be reformulated to allow in- 
dependent specification of Un/ut and ia/^t* Additionally, the new formula- 
tion had to be shown reasonably accurate insofar as ability to predict rotor/ 
turbulence noise is concerned, so that reasonable confidence could be estab- 
lished in the parametric study trends and results of the NASA-Ames 40x80 
wind tunnel evaluation. The following section describes the formulation of an 
anisotropic turbulence model for rotor/turbulence interaction predictions. 
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4.0 ANISOTROPIC TURBULENCE SPECTRUM ANALYSIS 


Recognizing that the assumption of isotropic turbulence is inadequate for 
modeling the turbulence spectra drawn into the fan duct, the next higher order 
in complexity is to assume that the turbulence is axisymmetric. Axisymmetric 
turbulence is characterized as having one preferred axis of symmetry, whereas 
isotropic turbulence has none. The axis of symmetry is assumed to be parallel 
to but not necessarily coincident with the ^main axial flow direction. The 
time-averaged velocity fluctuations (and length-scales) normal and parallel to 
the axis of symmetry may be different, whereas for isotropic turbulence they 
are equal. Adoption of an axisymmetric turbulence model, therefore, provides 
the flexibility to independently specify the axial and transverse turbulence 
intensities, as well as differing axial and transverse length scales. Be- 
cause only one axis of symmetry is assumed, the tangential and radial intensi- 
ties are also assumed to be equal, as are the tangential and radial length 
scales. 


The following development is based on the theory of Axisymmetric Turbu- 
lence given in References 10-12. The axial direction (subscript 1) is 
assumed to be the axis of symmetry. An effort has been made to develop the 
results so that they will reduce to the previously assumed case of isotropic 
turbulence when the proper limits are taken. The rotor/turbulence interaction 
model treats the rotor as a two-dimensional cascade in the axial-tangential 
plane, from which it follows that only the axial and tangential turbulence 
quantities are required. Denoting the tangential direction by subscript 2, 
the turbulence spectra ♦i2» ^22 required as functions of wave- 

number components k^, k2, and k3. The corresponding velocity correla- 
tion functions are R^x, Ri2» and R22» which are functions of spatial sepa- 
ration distances xx (axial), X2 (tangential), and X3 (radial). 

Chandrasekhar (Reference 11) showed that the general form of a symmetric, 
second-order, axisymmetric, solenoidal (divergence- free) velocity correlation 
tensor is as follows: 





( 1 ) 
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where 


’ij ■ \ [Hjk"i * •lu * “j’a)] 



and where Qi and Q 2 are arbitrary functions of |xij and xt ■ 1 ^x 2 ^ + X 3 ^. The 
symbol is the usual alternating tensor and is the Kronecker delta 

tensor. The above result is a generalization of the isotropic result 


■ lli 


R. . ■ x.x . + g 4 . . 

IJ I j ® *IJ 


where f and g are functions of r ■ /^l^ * 3 ^ 


( 2 a) 


and 


g 


2 r 


d / 2-v 

f) 


( 2 b) 


Previously (References 2-4), the isotropic result had been used, with f 
assumed to be of the form u^ exp (-r/jt). 


From Equations 1-3, expressions for the correlation tensors R]^^, R]^2* 
and R 22 can be derived, as follows: 

9 Q 1 

*^11 - -2Ql ' *t <3a) 



jQl 



21 “ *2 




2 Qi -Q 2 

aQl 

- 2 x, — - 2 

1 3 Xj 

3Ql 

3Q2 
“ *3 9x3 

2 Xj^Xj 

3^Qi 2 

-- X 

2 a'^l 

V : ■■ 

3 * 13 X 3 1 

2 

3 X 3 

*3 

3*i2 


(3b) 


(3c) 
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Whereas the isotropic turbulence correlation tensors, Equation 2, are charac'- 
terized by one function f(r), the axisymmetric turbulence correlation tensors, 
Equations I and 3, are characterized by two arbitrary functions Q]^ and Q 2 
of variables |xi| and Xf 

The turbulence velocity spectra corresponding to the velocity correlation 
functions are given by the Fourier transform of the correlations, as follows: 


m 

♦ j'ik) ■ I Rji(x) exp (-j k»x) dx 

ij J ^ - - - 


(4) 


Thus, if suitable functions Qj^ and Q 2 can be found and substituted into Equa- 
tion 3, relation 4 can be applied to the result to arrive at velocity spectrum 
functions 411 , ♦IZ. ♦ZZ* 

Certain measurable correlation functions can be derived as limiting forms 
of Equation 3. For example, the axial cross correlation of the axial turbu- 
lence velocity is given by 


Rjl (Xj,0,0) - -ZQj^ (xj,0) 


(5a) 


Also, the transverse (tangential) cross correlation of the tangential turbu- 
lence velocity is given by 


R22(0.X2.0) * -2Qj^(0,X2) -Q2(0»*2^ 


(5b) 


It is common practice to curve-fit the measured cross-correlations of Equa- 
tions 5a and 5b with functions of the type 


Rll(xjO,0) » u^^ exp (" 1 * 1 1 /fg) 


(6a) 


(6b) 


where Ug and U(; are the axial and transverse intensities, and tg and ft: 
the axial and transverse length scales, respectively. Equation 6a can also 
be written as an autocorrelation by assuming Taylor’s Hypothesis of frozen, 
convected turbulence. Thus, if the convection velocity in the axial direction 
is Ug, Equation 6a can be written as Ug2 exp(-Ug x/lg), where x is the auto- 
correlation time delay. 
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If expreisiont for and Q 2 found which satisfy Equations 5 snd 6 

in the proper limita, the above-described process for deriving fni ^ 1 . 2 , snd 
^22 carried out. Initially, forns for Q)^ and Q 2 vere assumed as 

follows: 

/ 

2 






* 't /*t ) ] 

(7 a) 




(7b) 


These expressions yield the empirical forms of Equation 6 for Rxi(xi,0,0) and 
R22(0 ,x 2,0) when substituted into Equation 3 and the appropriate limits are 
taken, per Equations 5a and 5b. 

Employing the assumptions for Q]^ and Q 2 given by Equations 7a and 7b, Che 
velocity correlation functions given by Equations 3a-c were derived, and Che 
corresponding turbulence spectra were derived by Fourier transform of the 
results, per Equation 4. The details of the derivations are omitted herein 
for brevity, but Che resulting spectrum expressions are as follows: 


.iVeV /k 2, 2 2 


♦12 


♦22 


2 2 
a t a 

w A X 


a 


(kita)(k2lj.) 


^ * f ] 

* A ^ f 

2X A / 


(8a) 


(8b) 


(8c) 


where 


A‘ 

X 


2 2 2 2 2 2 

‘ 1 + k/t/ * kj2,^2 ^ 

(9a) 

«,/t. 

(9b) 

“t^“a 

(9 c) 
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LboMlIbte.. .. 


The rotor/turbulence noite prediction requires the integral o£ the 
spectra over the spanvise wave nuaiber k 3 , defined as follows: 



(ki.ka) 


M 



^ • (kj » k2» kj) dk^ 


( 10 ) 


Utilizing Equations 8 and 9, application of Equation 10 yields the following 
relations for fn, ?22* 


■ •^^^(’‘'2 ‘t ** o ) 


(lla) 


3ij^^u ^ r 1 

♦l2 ^ 5 [^*''l^a^^‘‘2*t J 


(lib) 


[3(k^t^)^2-X^) 
4s L 


(lie) 


where 


. aJ (2p2 


2 2 2 2 
1 . k^^/ t . 


2+X^) 


] 


When expressions ll were progranuned and installed in the rotor /turbulence 
noise computer program, problems were encountered when attempting to compute 
cases where the axial-to>tangential turbulence length scale ratio X was 
greater than ~3. The source of this difficulty was traced to negative values 
of |22» which were physically unrealistic. Examination of the expression 
for I 22 gi-ven in Equation 11 shows that, in order to avoid negative values, 
we must have 

3(kjl^)^2-X^) + A^(2p^-2+X^') > 0 

17 




Now for the range 05(ki**)5 - , the ratio varies over the range 

from 0 to If At the lower limit» the above inequality gives 

1 \? « 2 + > 0 

which holds for all X «« long as V>1, At the upper limit, the above inequal'* 
ity gives 

2 2 
-f 2 - X > 0 

which only holds as long as 
X^ < 2 + 

For a typical velocity ratio M ■ 2, the maximum allowable value of length 
scale ratio is thus X « 2,45. Thus, for X > 2,45, negative values of ?22 
should occur, explaining the computational difficulties which were encountered. 

The source of this difficulty was found to be the assumed forms for Q], 
and Q 2 given by Equation 7, An additional constraint on the selection of Qj^ 
and Q 2 must be enforced which insures non-negative autospectra and 422. 

An extensive analysis of this problem was carried out to determine the 
forms of the spectrum functions which are non-negative for all practical 
values of X and V. It is shown that (see Reference 12) the spectra can be 
expressed in the form 

4^1 (k) - k 3 ^) F (k) (12a) 

*12^-^ * - kj^k2 F (k) (12b) 

422(k) - (kj^ + kj^) F (k) + G (k) (I2c) 

where F(k) and G(k) are arbitrary functions of ki^ and kj;2 ■ k 2 ^ + k 3 ^; in 
order for 4]^]^ and 422 to be non-negative, they must satisfy the following 
constraints : 

F (k) >0 (I3a) 

F(k) + (k2 + kj) k"^ G(k) >0 (I3b) 
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Note that # 1.1 and #i2 given by Equation 8 can be put in the fonn given 
by Equation 12 and additionally oieet the requireawnt F({c)>0. The exprea- 
sion for ^22 given by Equation 8 can also be put in the form given by Equa- 
tion 12, but doea not meet the conatrainta of Equation 13. Particular func- 
tional forma for F(k) and G(k) have been derived which aatiafy inequalities 
13 , and they are aa followa: 


F(k) - 


2t u ^ 

a t a 

TT5 

» A 


[k) 2 


G(k) - F(k) I 2P^ - I - l/X^ 


] 


(14a) 


(14b) 


The above expreaaion for F(k) is identical to that implied for Equation 8, 
but the relation above for G(k) la much aimpler. If Equationa 14 are sub- 
stituted itsno Equation 12, it is easily derived that 422 remains non-negative 
as long aa 


2 2 
x2> Jl. or ^ 


(15) 


Since, for the rotor/turbulence interaction problem, u^^ < u^2 and > i^, the 
above constraint is usually always met. 

To obtain the two-dimensional spectra, expressions 12 and 14 are substi- 
tuted into 10, Integrating over spanwise wavenumber, to obtain the following: 


♦ (k k ) ■ ^ ^ (3 k ^ + A 

♦ jlUj.Kz; . . 5 *^2 t ^o ^ 

4v A 


31 ^u ^ 

Air A 




* 

a t a 
4v A ^ 


[’“i 


(i/x^) 


+ (2P^ - H\h A ^ 
o 


(16a) 


(16b) 


(16c) 
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It can be seen that these expressions are the same as those previously de- 
rived, Equations 11, with the exception that a term (2-v^) $22 Equa- 

tion 11 is now replaced by (l/x^)> The requirement for $22 
negative is the same as that given by Equation 1S> 

It may be noted in passing that the final form for G(k) given herein 
(Equation 14) is very similar to one expression suggested by Sreenivasan, 
Reference 12, i.e., 

G(k) - F(k^ - 2] . 

Expressions 16 were programmed and installed in the rotor/turbulence 
interaction program, and diagnostic calculations of narrowband (1 Hz band- 
width) PWL were carried out for the same case (Rotor 11 at 54Z speed, from 
Reference 6) as was shown in Figure 3. Axial and transverse length scales 
were parametrically varied, and results of these calculations are shown in 
Figures 4 through 8. Figure 4 shows the portion of the spectrum centered 
around blade-passing frequency (6PF), for a fixed transverse length scale of 
one blade spacing. Figure 5 shows the corresponding results at 2 x BPF. 

Note that the spectra are very peaky for axial scales t^/s > 50, and that 
the second-harmonic spectra are broader than the fundamental spectra. The 
fall off in peak tone level with harmonic number (n) is shown in Figure 6. 

The curves have similar shapes for ig/s > 10. The increase in tone peak level 
at BPF with tg/s is shown in Figure 7. Above about i^/b > 50-70, the tone peak 
increases linearly with Jlg/s. The effect of transverse length scale varia- 
tions, wtih axial scale held fixed, is shown in Figure 8. Smaller it/s in- 
creases the peak level and makes the spectrum narrower, but the peak appears 
to maximize for this case at about t^/s * 0.1. 

From the results shown in Figures 4 through 8, it can be expected that, 
for typical values of tg/s s; 50-100 and i^/a z 0. 5-1.0, the rotor turbu- 
lence interaction will produce spectra which are very peaky around blade- 
passing frequency and its harmonics, and will have broadband levels which are 
substantially below the peaks, by as much as 40-50 dfi or more. A complete 
spectrum for one combination of ta/s and i^/s typical of that measured in 
Reference 6 is shown in Figure 9. This result confirms that the expectation 
is correct. 
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It was observed that the spectrum peaks at blade-passing frequency be- 
came higher but narrower as i^/s increased, Figure 4. It was postulated that 
the total tone PWL (integral over a prescribed bandwidth) might become con- 
stant above a certain value of A^/s and remain so with any further increases 
in 4^/s. Numerical Integrations were performed over a 100 Hz bandwidth 
atound BPF and its harmonics for the results shown in Figures 4 through 8> 

The total PWL in a 100 Hz bandwidth does indeed level off with increasing 
i|i/s, as the results in Figures 10 and 11 show. 

Because the spectra can be so peaked for large length scales (see Fig- 
ure 9), it is no longer sufficient to compute 1/3-octave levels by evaluating 
the center frequency level and multiplying by the appropriate bandwidth, as 
had been done in the past (e.g.. Figure 3). The 1/3-octave levels must be 
evaluated by integrating the narrowband spectrum over the bandwidth. As Fig- 
ure 9 shows, this involves computing many points within a given frequency band 
in order to obtain an accurate representation of the spectrum shape, since the 
spectrum drops off so rapidly on either side of the BPF harmonics. This pro- 
cedure can be costly and time-consuming, and so alternate methods for evaluat- 
ing 1/ 3-octave levels accurately were investigated. 

A method was finally developed whereby the 1/3-octave levels are computed 
by evaluating the narrowband levels at selected points in the 1/3-octave band. 
The integration is performed by curve-fitting the computed points with simple 
mathematical functions whose integrals are available in closed form. This 
procedure is illustrated in Figure 12, where a BPF harmonic exists at some 
arbitrary location within the 1/3-octave band. Narrowband levels are computed 
at the peak and at both ends of the band, as shown in Figure 12a. 

Considering one-half of the spectrum, e.g., f > nxBPF, as shown in Fig- 
ure 12b, the peak level is designated PWLq and sideband value as PWL^. These 
hcye corresponding power values of Pq and P]^, respectively. The sideband 
frequency is denoted by f]^. From the shapes of the curves shown in Figures 4 
through 9, a logical choice for a curve shape is 
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Figure 10. Predicted Effect of Axial Length Scale on Integrated 100 Hz 
Bandwidth Tone PWL. 





Narrowband PWL, dB 
















This can be Integrated over the range Af ■ (f i-fg) to give 


**t " iniPp/Pj) 

vhere P^ is the total power in the interval Another logical choice 

for curve-fitting the spectra is 

P - P^/[l ♦ a^(f-f^)2] 

Again integrating over the interval Af “ ^l^fo and evaluating the constant 
a from the sideband values of P and f-fg, we obtain 

These two expressions were compared with numerical evaluations over 100 Hz 
and 1/3-octave bands, and results showed that either of the above appruxima^ 
tions 17 or 18 gave accurate results for aPWL ■ PWho-PWL], <20-25 dB. For 
APWL > 20, it appears as though the average of 17 and 18 gives the best com- 
parison with the precise numerical integrations. This is shown in Figure 13, 
where the decibel equivalents of 17 and 18 are plotted and compared with 
numerical integraton points. The average curve of Figure 13 is therefore 
used in the present computer code for evaluating 1/3-octave spectra. 

The final formulations discussed above permit evaluation of rotor/ 
turbulence interaction noise when the turbulence properties are known at the 
fan face (rotor leading edge plane). The turbulence properties Ug, u^, ig, 
and At can be specified independently, consistent with measured values. How- 
ever, sometimes the properties may only be known upstream of the fan face at 
a location where the flow still undergoes a contraction prior to entering the 
rotor. A method for estimating the turbulence spectra for axisymmetric tur- 
bulence passing through a contraction is therefore required. The sudden con- 
traction theory of Ribner and Tucker, Reference 5, is utilized herein. 


31 






Let Che posCconCraccion speccra and wave nuaber be given by end 

Che pteconCracCion apecCra and wave nuaber be given by tij(h)< C denoce 

Che conCrecCion racio, defined as Che racio of downscresaoCo-upaCceam flow 
velocicy, and define c ■ L/C^. Given k ■ (k|, K2, K3), Che corresponding 
preconc race ion wave ntaiber is given by k ■ (ki, k2i k3> ■ (Ck^, K2//?» 

*3/1^). An axisymaeCric conCracCion is assuaed, as discussed in Reference 5 . 

Per Reference 5 , Che posC conCracCion specCra are given by Che following 
formulae, in Ceras of Che preconCracCion specCra 



24ii(k) kj^^d-e) 

a j j 

c k/ + kj" ♦ kj" 


(ekj2 + kj^ + ^ 


( 19 ) 


♦l2<S> 


7S * 


[♦ll(k) kj ♦ ♦i2(k) k^] k^(l-c) 

V.k/.k/ 


♦ll(k) k^^ kj (l-e)^ 

(ck ^ + k ^ + k3^)^ 


( 20 ) 


♦22<!S> 




29i 2^-^ kj^k2(l-e) 
2 2 2 
^ h * V 


♦u<!s> 

n 9 2 2 

(ek^^ ♦ k^** + k^n^ 


( 21 ) 
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The corregponding Cwo-dimenalonal specCre, defined as the integral over 
K 3 of the above three-dimensional spectra, are then given by 


m 

*i.(ki,k2) • 2 /C j 


OCj.'Cj.'Cj) dkj 


(22) 


where ^|j_j(k) are assumed to be even functions of k 3 and the relation dtc 3 • 

4 k 3 was used. The above Integrations (Equation 22) are carried out numer- 
ically. Let 


i. k- ■ A tan 6 
t ^ o 


% ‘ a *^a t '^2 


(23a) 


(23b) 


then it**k 3 ■ Aosec^OdO, and the factor A“^ in the expressions ♦ij(k), Equa- 
tions 12 and 14, can be written as 


.-6 


(1 * , ,^2,^2 j 

(“o' * *tS') 


-3 


A -6 -6- 

A sec 6 
o 


so that the integrals over k 3 can be written in the form 
-If/ 2 


2/c f 

lA^ J 

to n 


[ 


'J cos^ 0 d 6 


The interval O<0< ir/2 is divided into 18 equal increments (5*), and 
Simpson's Rule is used to evaluate the interval. Equations 12 and 14 can be 
written in the more general form 


'<’ij (k) 


• i 2 2 

« ^ ^ u ^ 

1 - „ ■ ? - r., cos^ e 

"J 


(24) 
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where 


2 2 2 2 
,, - k„ i + k, r 

(25a) 

11 2 t 3 t 


12 “ " (kiA3)<k2«-t)/^ 

(25b) 

22 “ (2|i^-l/X^) 

(25c) 


Equations 19 through 25 provide all the necessary relations needed to compute 
the postcontraction, axisytimetric turbulence spectra. 


5.0 SURVEY OF INLET TURBULENCE PROPERTIES 


A literature survey of published data on fan inlet turbulence properties 
was carried out. A sumnary chart of the data reviewed is given in Table 1. 

The table lists the reference^ vehicle, environment (test stand, chamber, wind 
tunnel, etc.), the approximate values (or range of values) of turbulence 
properties measured, and some of the important fan geometric parameters. From 
a review of the references listed in Table 1, an "expected ensemble average" 
set of turbulence properties was deduced for a typical anechoic chamber, a 
typical outdoor test stand, and the NASA^-Ames 40x80 wind tunnel. Expected 
average values for in*-*flight conditions were derived from the results of 
flight test correlations of atmospheric turbulence summarized by Houbolt in 
Reference 25. 

ii table of typical expected values of inlet turbulence velocities and 
scales are given in Table 2. The chamber values rely heavily on the infor- 
mation in Reference 6, since all four quantities were measured, as well as 
the fan noise characteristics. Additionally, some of the earlier data was 
associated with peculiar facility configurations and/or had inadequate in- 
strumentation and data acquisition equipment for measuring long length scales 
(see Reference 6). The outdoor test stand expected values in Table 2 are 
weighted heavily by the data of Hanson (References 1 and 16), as well as pre- 
viously unpublished taken on a fan engine at the General Electric 

Company Peebles, Ohio, test facility. 

For the flight case, no fan inlet measurements of turbulence were avail- 
able. However, the atmospheric turbulence properties can be used (Reference 
25) along with contraction ratio corrections as outlined in the previous 
section. The contraction ratio is assumed to be equal to the fan face axial 
velocity divided by flight velocity, i.e., C = ^a^^o* 

The NASA-Ames 40x80 wind tunnel properties are taken from Reference 20. 
The measurements were made 20 feet upstream of the fan inlet, and no fan face 
measurements were available. In order to utilize the upstream measurements, 
several steps must be taken. First, both tangential and axial turbulence 
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Table 1. Summary of Fan Inlet Published Data 
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Table 2. Summary of Average Expected Inlet Turbulence 

Properties for Turbofan Engine Test Facilities. 




velocities are given, but only axial length scale was measured. Hence, 
the transverse length scale must be deduced from the ratio of tangential 
to axial turbulence velocity and the sudden contraction theory of Reference 
5. This amounts to assuming that somewhere upstream of the measurement point 
the turbulence was isotropic, and that the eddy elongation and cross section 
contraction correspond to the measured velocity ratio. Second, the turbulence 
at the measuring point is assumed to undergo another "sudden" contraction from 
the measuring point to the fan face, and the relations developed in the pre- 
vious section are used to compute the fan face spectra from the measuring 
point spectra. The contraction ratio in this case is assumed to be the ratio 
of fan face velocity to tunnel velocity* i.e., C <■ Uq/Vq. 


6.0 DATA/THEORY COMPARISONS 


Having developed an axisymmctric turbulence spectrum model, extensive 
data/thcory comparisons were carried out to verify that the new prediction 
model gave reasonably accurate predictions of rotor/turbulence interaction 
noise. The data of Reference 6 were first compared, since it contained 
measurements of both turbulence properties and farfield acoustic spectra. 

The fan stage, a 0.508 m (20 in.) diameter fan with A4 rotor blades and 86 
stator vanes, was tested with and without a honeycomb /screen turbulence control 
structure around the fan inlet. The purpose of the control structure was to 
suppress the incoming turbulence and eliminate the rotor/turbulence noise. 

For the purpose of this study, however, the control structure produced a 
different type of turbulence at the fan face and hence a different rotor/ 
turbulence interaction. It is of interest to sec if the prediction model cor- 
rectly predicts the effects of the control structure. 

The noise predictions were made using the rotor blade geometric proper- 
ties and mean flow conditions at the root-mean-square radius of the inlet an- 
nulus, i.e., at a "pitchline** radius given by 

rp = /(r§ + rg)/2 

The lift coefficient (used to determine the quadrupole strength) was evaluated 
based on incidence angle at the pitchline rp, rather than from pressure 
ratio (Reference 3), since it was felt that leading edge loading was the primary 
influence on forward-radiated noise. The quadrupole portion of the prediction 
model has a singularity at M,. 1, where is the rotor inlet relative Mach 

number , 

My » + M^ 

Therefore, predictions had to be limited to those speeds for which My < 1. 

The design tip speed for this fan is Uj- = 427 m/s (1400 fps), so comparisons 
were limited to speeds less than 86% of design. Actually, the quadrupole pre- 
diction model is based on linear flow theory which is not valid in the tran- 
sonic regime, so accurate results should not be expected for My > 0.9, which 
corresponds to 78-80% speed. 
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Computed and measured 1/3-octave spectra are shown in Figures 14 through 
18 for the HASA Rotor 11 fan stage described above, without the turbulence 
control structure (TCS) in place. Corresponding spectra with the control 
structure in place are shown in Figures 19 through 22. A summary of the data/ 
theory comparisons made is listed in Table 3. Only inlet quadrant PWL spectra 
arc shown. 

Considering first Figure 14a, which is the same case shown in Figure 3 
(using the isotropic turbulence model), it can be seen that the prediction of 
blade-passing frequency tones and harmonics is in good agreement with the data 
for the new axisyianetric turbulence formulation. The results are a tremendous 
improvement over the isotropic predictions of Figure 3. Note also that the 
broadband rotor/turbulence levels are substantially below the measured levels, 
as would be expected from the narrowband prediction shown in Figure 9. 

The blade-passing frequency for the case shown in Figure 14a is app’-oxi- 
mately 6300 Hz, and a BPF harmonic is predicted to occur in every 1/3-octave 
band above 20,000 Hz. The predicted levels above 20,000 Hz are in good agree- 
ment with the measured levels, although the predictions do not contain air 
attenuation effects. This would tend to lower the predictions at the high 
frequencies, by as much as 11 dB at 80 KHz. Examination of narrowband measured 
spectra (e.g.. Figure 48a of Reference 6) shows that the BPF harmonic tones 
dominate their 1/3-octave band up to the 5th harmonic. Beyond n*5, the broad- 
band noise dominates the 1/3-octave band, and this is consistent with the 
predicted curve of Figure 14a if air attenuation effects were to be accounted 
for. Similar remarks apply to the other speeds and discharge valve (DV) set- 
tings shown in Figures 14 through 18. 

It is observed that the theory overpredicts the noise at high frequencies 
at 742 speed, where Che relative Mach number is rather high, ~ 0.85 at the 
pitchline. At 80% speed, Che overprediction worsens, and M^ >_ 0.9 for this 
speed . 

Considering now the cases with the TCS in place, Figures 19 through 22, a 
comparison of the meas'un‘d average fan face turbulence properties with and 
without the TCS is shown in Table 4. These values represent approximate 
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Figure 19. Comparison of Predicted and Measured Inlet PWL Spectra for NASA Rotor 11 
Fan Stage with TCS at 54% Speed. 
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Figure 21. Comparison of Predicted and Measured Inlet PWL Spectra for NASA Rotor 11 
Fan Stage with TCS at 74% Speed. 
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averages of all the data taken in Reference 6. A range of length scale ratios 
is shown in Table 4 with the TCS in place, because the tangential length scale 
was found to be smaller than was measurable with the instrumentation employed 
by the authors of Reference 6. The low limit corresponds to assuming that &|- 
is the same with the TCS as without, while the high limit corresponds to as- 
suming that the scale ratio X ■ la/^t same with TCS as without. The 

actual scale is probably somewhere in between the two limits. 

The data/ theory comparisons for the cases with TCS are shown in Figures 
19 through 22. The shaded band for the predictions corresponds to the range 
of length scales assumed, 90 £ ta^^t £ general, the case ia/^t * 

corresponds to the high limit of the shaded band for frequencies below 3 x BPF, 
and the case ta/^t * corresponds to the high limit for frequencies above 
3 X BPF. The general observation to be made about the results shown in Figures 
19 through 22 is that the "predicted noise reduction due to addition of the TCS 
is at least as much as was measured, if not more. The predicted spectra indi- 
cate that residual rotor /turbulence interaction noise still exists at BPF for 
the DV * 0 (open throttle) cases, whereas the DV = 1.27 (closed throttle) cases 
do not appear to have any appreciable rotor/turbulence noise contribution. 

A comparison of measured and predicted BPF tone levels, i.e., the 1/3- 
octave band level containing the tone, is shown in Figure 23. The underpre- 
diction at 60% speed without TCS (see Figure 15) is due to the BPF being right 
at the edge of the band. The agreement is good except at high speeds where 
transonic conditions exist and the theory is not really expected to be accurate. 

Data/ theory comparisons were also carried out for a scale model variable- 
pitch fan tested in the same anechoic chamber (General Electric Research and 
Development Center, Schenectady, New York) as was the NASA Rotor 11 fan stage 
reported in Reference 6. Inlet arc 1/3-octave PWL spectrum comparisons are 
shown in Figures 24a-f, the data being taken from Reference 26. This fan has 
18 blades and 33 vanes, and has a design tip speed of 306 m/ s (1000 fps). It 
is a variable-pitch design and consequently has a radially constant solidity 
(chord/spacing), as opposed to conventional fixed blade fans which usually 
are designed to have radially constant chord and variable solidity. It also 
has significantly fewer rotor blades than the NASA Rotor 11 fan stage (18 
versus 44), and has a subsonic design tip speed. 
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Table 4. Summary of Assumed Turbulence Quantities 
for GE Schenectady Anecholc Chamber Fan 
Noise Facility. 


Table 5. QCSEE Scale Model Variable Pitch Fan 
Data/Theory Comparison Summary Chart 


Case No. 


% Np 

DV 

TCS 

(Mr)tip 

60.2 

7.75 

No 

0.620 

70.4 

7.75 

No 

0.731 

80.4 

7.75 

No 


90.5 

7.75 

No 


95.5 

7.75 

No 

1.063 

99.2 

7.75 

No 

1.074 


<Mr)pitch Figure 


0.514 

0.607 

0.706 

0.811 

0.866 

0.903 




































The predictions shovm in Figures 24a**£ were done in the same fashion as 
was done for the NASA Rotor 11 fan> i>e,, section properties at the rms pitch- 
line radius were used* ' The turbulence properties of Table 4 were used as in- 
put, since the same facility was used. The data was taken without a TCS, and 
only along the QCSEE engine sea level operating line. The data shown in 
Figures 24a-f were taken with a hardwall, low Mach number inlet installed. 
Again, input lift coefficients for the predictions were based on pitchline 
incidence angles rather than rotor pressure ratio. A summary of comparison 
points is given in Table S. 

The same observations made previously for the Rotor 11 data/ theory com- 
parisons apply for the QCSEE Variable Pitch fan also. The EPF tones are pre- 
dicted quite well, the rotor/ turbulence broadband noise contributes nothing 
to the observed broadband levels, and the theory overpredicts at transonic 
relative Mach numbers. A summary of the QCSEE fan stage data/ theory compari- 
sons for BPF harmonics is shown in Figure 25. 

It can be seen from Figures 24 and 25 that the predictions indicate only 
a small contribution of rotor/turbulence noise at 3 x BPF and higher, espe- 
cially at the lower speeds. The progresively larger contributions at higher 
speeds is consistent with the change in measured spectrum shape with speed, 
i.e., the BPF higher harmonic peaks progressively stand out more above the 
broadband levels as speed increases. It is possible that the input lift coef- 
ficients based on incidence angle are too low, and this would yield lower pre- 
dicted levels at frequencies of 2 x BPF and higher. The few narrowband samples 
given in Reference 25 do show, however, that the tones do not contribute 
significantly to the 1/3-octave levels above the third harmonic, supporting 
the predicted trends of Figure 24. 

An attempt was made to predict the blade passing frequency tones for the 
Pratt & Whitney JT15D engine, and compare the predictions with the experimen- 
tal results recently obtained* in the NASA-Ames 40x80 ft wind tunnel. Tur- 
bulence measurements were made previously, reported in Reference 20, approxi- 
mately 20 ft upstream of the engine inlet, but no data was available on the 


^Contract NAS2-8675 in progress, unpublished data. 
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\ X Bl>F 



Figure 25. Comparison of Measured and Predicted a x BPF 
Tone PWL for QCSEE Scale Model Fan Stage. 










turbulence at the fan face. Linear sudden contraction theory. Reference 5, 
was therefore enployed to infer the turbulence properties at the fan face, 
based on the reported ncasurements upstreasi of the inlet. The measurements 
in Reference 20 showed the axial turbulence intensity in the tunnel test 
section to be '*0.2Z, independent of tunnel speed. Likewise, the trans** 
verse turbulence intensity was found to be '^0.55%, also independent of 
tunnel speed. 

Using an engine test condition corresponding to a fan tip speed of 1020 
fps (the lowest tip speed for which data was taken in this series’*'), the fan 
face velocity is approximately 394 fps. For tunnel speeds of 11 knots and 80 
knots, this implies a turbulence contraction ratio (or speed ratio) of 21 and 
3, respectively. Applying the theory of Reference 5, the turbulence properties 
computed at the fan face are as follows: 

Turbulence Component Vq ■ 11 kts. ■ 80 kts. 


Axial intensity 

0.00084X 

0.032X 

Transverse intensity 

0.12X 

0.32X 

Axial length scale 

7.56 ft 

3.97 ft 

Transverse length scale 

0.015 ft 

0.15 ft 


It can be seen from the above tabulation that the estimated turbulence inten- 
sities are extemely small; the corresponding blade passing tone PWL levels 
were predicted to be 99.0 dB and 107.3 dB, for 11 knots and 80 knots, respec- 
tively. The corresponding measured values were 123.8 dB and 117.5 dB, for II 
knots and 80 knots, respectively. 

The above results would seem to imply that rotor/turbulence noise in the 
40x80 ft wind tunnel is not a significant contribution. Howevet, the fact 
that the measured noise decreases by 6.3 dB from a wind velocity of II knots to 
80 knots indicates some sort of flight "clean-up" effect is taking place. In 
fact, examination of 20 Hz narrowband spectra of the far field microphone 
data shows that the II knot wind data exhibits skirted peaks at blade passing 
frequency typical of a turbulence-modulated tone. The 80 knot data, however, 
shows much narrower peaks . 

*Ibid. 
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One possible explanation for the above inconsistencies is that the 
sudden contraction theory may not be an adequate representation of what 
happens to the turbulence in a wind tunnel as it contracts in passing through 
the engine inlet. There is some evidence, reported in Reference 24, that the 
transverse intensity does not change at all from upstream to the fan face. 
Further, there is also some evidence* vht?t axial intensity does not change 
appreciably from the tunnel test section to the fan face. 

Examination of the experimental measurements of wind tunnel contraction 
on free stream turbulence made by Uboroi, Reference 27, shows that the reduc- 
tion in axial turbulence intensity across a contraction is not nearly as great 
as the theory of Reference 5 would predict for large contractions. The 
corresponding change in transverse turbulence intensity is overestimated by the 
theory, but the error is not nearly as great as that for the axial intensity. 

A summary of these trends are shown in Figure 26. 

The differences shown in Figure 26 would affect the 11 knot predictions 
but would have no significant impact on the 80 knot predictions, since the 
contraction ratio Uq/Vq (fan axial velocity/tunnel velocity) is only about 
3:1 for the 80 knot case. The 11 knot condition, however, was experimentally 
obtained by turning off the tunnel fans and allowing the engine to pump the 
tunnel naturally. This condition could therefore produce turbulence charac- 
teristics which arc significantly different from those measured in Reference 
20 with the tunnel fans in operation. Extrapolating the turbulence data of 
Reference 20 to 11 knots (the minimum tunnel speed for which data were taken 
was -40 knots) could therefore yield an incorrect assessment of the turbu- 
lence properties for the 11 knot case, contributing to the poor agreement 
between experiment and prediction. 

Another source of the differences obtained between predicted and measured 
BPF tone levels is that the predictions were originally made utilizing only 
the dipole source contribution, and the quadrupole source contribution was 
omitted. This was done because there was insufficient information available to 


^Personal communication, C. Feiler and L. Shaw of NASA Lewis Research Center. 
Results obtained in Lewis 9x15 wind tunnel. 








GE on the aerodynamic characteristics of the JT15D fan rotor. If, however, 
it is assumed that its loading characteristics are similar to current high 
speed fan designs in the GE family of engines, a reasonable estimate of the 
blade loading (or lift) coefficient can be made, and this is required to pre- 
dict the quadrupole contribution to the BPF tone level. Assuming a typical 
lift coefficient of * 0.45, the quadrupole contribution at Vq * 80 knots 
was found to be equal to the dipole contribution, thus increasing the pre- 
dicted level to 110.3 dB, within "7 dB of the measured level of 117.5 dB. 

Another prediction was also made for the 80 knot case where the length 
scales were assumed to obey the sudden-contraction theory behavior as far as 
elongation of the axial length scale and contraction of the transverse length 
scale is concerned. The turbulence intensities, however, were assumed to be 
the same at the fan face as was measured upstream, as was observed by Hodder 
in Reference 24. This prediction yielded a BPF tone level of 116.0 dB, within 
1.5 dB of the measured value of 117.5 dB. It is recognized that Hodder *s re- 
sults may not apply to the 40x80 wind tunncl/JT15D engine application, since 
they were obtained in a much smaller tunnel on a low speed fan. Nevertheless, 
the majority of the evidence collected thus far indicates that the sudden- 
contraction theory is inadequate for estimating fan face turbulence Inten- 
sities from measurements made far upstream of the inlet. The empirical ad- 
justments of the sort described above seem warranted and do yield reasonable 
agreement between measured and predicted noise levels. 

To test the ability of the prediction model to evaluate wind tunnel tur- 
bulence noise for a fan, predictions were made of the BPF tone PWL (inlet arc) 
for the 15-blade fan tested in the NASA-Lewis 9x15 wind tunnel, the results of 
which are given in Reference 21. In this case, turbulence measurements were 
made inside the fan duct close to the rotor, so the question of the inadequacy 
of the sudden-contraction theory for predicting turbulence intensity changes 
is not an issue. For the static (Vq^O) case, the measured transverse and 
axial intensities were used to estimate the effective contraction ratio, using 
the sudden-contraction theory, and this contraction ratio was used to estimate 
the transverse length scale only , based on the measured axial length scale. 

For the wlnd-on case, Vq = 80 knots, the actual contraction Uq/Vq was used to 
estimate the transverse length scale. 
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The resulting predictions of BPF tone PWL versus fan speed for both 
static and wind-on conditions is shown in Figure 27 along with the measured 
values. It can be seen that the agreement is, on the average, reasonably 
good, considering the approximations and assumptions made, lending further 
support to the validity of the basic rotor/turbulence interaction predic- 
tion model . 

Based on the analysis and review of existing data described above, it is 
concluded that the current rotor/ turbulence interaction noise model is adequate 
for predicting wind-tunnel turbulence/rotor interaction noise when the fan face 
turbulence properties are known. When only upstream turbulence properties are 
given, the sudden-contraction portion of the turbulence prediction may be in- 
adequate for large contraction ratios. For tunnel velocities of 80 kts. or 
less, this could imply an underprediction of the BPF tone PWL by about 7-8 dB, 
if the prediction model is used as-is. With the use of an empirical correction 
to maintain constant fan face turbulence intensity independent of tunnel speed 
or fan speed, the prediction will be fairly accurate. Whenever contraction 
ratio is close to unity, i.e., V^^Ua, either method of prediction will 
yield adequate results. 
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7.0 PARAMETRIC STUDY 


Having established that the present rotor/turbulencc interaction theory 
adequately predicts rotor/turbulence noise (for subsonic relative Mach numbers) 
over a wide range of fan operating conditions, fan geometries, and inlet tur- 
bulence conditions, a parametric study of rotor turbulence noise characteris- 
tics was carried out. It was shown in the previous section that rotor/tur- 
bulence noise is characterized by strong, narrow peaks at blade passing fre- 
quency (BPF) and its harmonics, and contributes an insignificant amount of 
broadband noise to the total fan noise spectrum. Therefore in the present 
parametric study, predictions are made only at blade-passing frequency and its 
harmonics. 

When only the BPF tones and harmonics are of interest, a considerable 
simplification of the prediction equations can be made when Ig^/S » 1. This 
simplification is described In detail In Appendix A, and permits calculation 
of the total sound power in the vicinity of BPF, and for each harmonic thereof. 
This method was found to be quite accurate when compared to l/3-octave cal- 
culations, over a wide range of fan operating conditons, inlet turbulence 
conditions, and fan geometries. The parametric study to be described In the 
following sections utilizes the simplified theory outlined In Appendix A. 

The prediction program requires as input the fan rotor inlet axial and 
rotational Mach numbers and rotor total pressure ratio. A composite plot of 
these fan parameters was made for several fan stages, including those discussed 
in the previous section. Figures 28 and 29 show plots of axial Mach number 
Mg and fan pressure ratio PR versus tip speed Mach number M^, respectively. 

It can be seen that many of the fans have a common flow versus speed (Mg versus 
M{^) characteristic, and a common work versus speed (PR versus M^.) charac- 
teristic. 

A common operating characteristic of Mg and PR versus M^. was therefore 
selected for the parametric study, and these are shown as the solid lines in 
Figures 28 and 29. This is henceforth referred to as the low-flow (LF) tran- 
sonic fan operating line. The QCSEE fan and NASA Rotor 55 fans, however, have 
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higher Mj| vecsu$ characteristics (Figure 28), aed this is related to 
their design tip speeds, idiich are subsonic. A limited number of calculations 
were therefore performed using the QCSEE flow characteristic (M^ versus M^), and 
this operating line is henceforth referred to as the high-flow (HF) subsonic fan 
characteristic . 

Several fan sizes were selected for study because it is known from theo- 
retical considerations that the turbulence scale-to-blade spacing ratio is im- 
portant. Hence different fan sizes operating in the same turbulence environ- 
ment may give different noise levels when the levels are corrected to a common 
fan inlet area. Fan diameters of 2.13 m (84 inches), 1.07 m (42 inches), and 
0.53 m (21 inches) were selected for study. They are henceforth referred to 
as full-scale (FS), half-scale (HS), and quarter-scale (QS), respectively. 

Blade number was also varied , values of Ng ” 38, 28, and 18 being selected 
for study, lliese values represent the range of fan types currently in use. 

For example the diameter/blade number combination 2.13 m/38 is close to the GE 
CF6 engine fan, and the combination 0.53 m/28 is the same as the P&W 41150 
engine fan. 

Blade tip solidity was held constant at (c/s)tip “1.3 for all LF tran- 
sonic fan cases and 1.0 for all HF subsonic fan cases, as it is usual design 
practice to set rotor tip solidity approximtely equal to design tip relative 
Mach number. A list of the fan geometry parametric variations investigated is 
shown in Table 6. 

Calculations of rotor/ turbulence noise were made for each of the configu- 
rations listed in Table 6, over the tip speed Mach number range 0.5 < < 1.0. 

The calculations were performed for several turbulence conditions, as follows: 

• Outdoor Test Stand (OTS) 

• AMES 40x80 Wind Tunnel, Vq « 80 kts (AWTl) 

• AMES 40x80 Wind Tunnel, Vq 180 kts (AWT2) 

• In-Flight, Vq - 180 kts (FLT) 
















Additional calculations were made for the quarter-scale fans in the Schenectady 
Anechoic Chamber (SAC). Noise levels were obtained at blade-passing harmonics 
n«l through 8 for all cases. Some parametric studies were also made of the 
effect of length scale and flight speed. 

A question was raised in the preceding section about the adequacy of 
using the sudden contraction theory (Reference 5) for estimating the fan face 
turbulence characteristics based on measured turbulence properties far upstream. 
In the present study, sudden contraction theory was used without alteration; 
however, estimates were made of the range or uncertainty band in predicted 
noise levels for the wind tunnel cases, and the uncertainty band was found to 
be less than 7 dB for the worst cases. A more detailed discussion of this 
topic is covered in the next section, NASA-AMES 40 x 80 WIND TUNNEL EVALUATION, 
but it is noted here that the 7 dB uncertainty band is not large at all com- 
pared to the uncertainties associated with turbulence variability, tone un- 
steadiness and day-to-day weather changes which affect the Outdoor Test Stand 
and Flight rotor/turbulence noise levels. 

7.1 COMPONENT SOURCE CONTRIBUTIONS 

The relative contributions of the dipole (unsteady lift) and quadrupole 
(steady loading/turbulence interaction) sources to the total predicted rotor/ 
turbulence noise was examined first. Figure 30 shows the contributions for 
the quarter-scale fans on the Low-Flow operating line, in the Schenectady 
Anechoic Chamber (SAC) environment. Note that the dipole source dominates the 
BPF tone for the 38-blade rotor, whereas the quadrupole source dominates for 
all but the lowest tip speeds (M^. > 0.6), for the 18-blade rotor. This is 
because the rotor blade aspect ratio is smaller for the 18-blade rotor, and 
the unsteady lift response function decreases as aspect ratio is reduced. Ihe 
lift coefficient, which determines the quadrupole level, is the same for the 
two blade numbers, since rotor solidity was held constant. For this parametric 
study, input lift coefficient was Cj^ ■ 0.45, based on Che low-flow operating 
line pressure ratio and flow versus speed. For this parametric study, using 
leading edge loading or incidence angle for computing Cj^ was not feasible 
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Figure 30. Component Source Contributions to BPF Tone PWL for 
Quarter-Scale Fans in Schenectady Anechoic Chamber. 







since the actual running incidence angle is an unknown function of aerodynamic 
design criteria and effectiveness of the design process, and accounting for 
this was beyond the scope of the present study. 

Figure 31 shows the component source contribution breakdown as a function 
of BPF harmonic number for the quarter-scale fans in the SAC turbulence envi- 
ronment. These results indicate that the higher harmonics of BPF are domini- 
nated by the quadrupole source, even for the high blade number Ng ■ 38. 

The fall-off with harmonic number n is very large for the dipole source compared 
to that for the quadrupole source. 

Similar trends are shown for the Outdoor Test Stand (OTS) turbulence en- 
vironment (Figure 32), the Ames 40 x 80 Hind Tunnel (AHT) environment (Figures 
33 and 34), and the Flight (FLT) environment (Figure 35). The above calcula- 
tions were repeated for the full-scale fan in the OTS, AWT and FLT turbulence 
environment, and these results are shown in Figures 36 through 39. Sample 
spectra, i.e., BPF tone PWL versus n, are given for both quarter- and full- 
scale fans, for each turbulence environment, in Figures 40 through 43. 

It can be seen from the results shown in Figures 30 through 39 that the 
trends of BPF tone PWL with Mach number are similar for SAC, OTS, and 
AWT, but the FLT condition yields a flatter curve. Also, the quadrupole 
contribution for the FLT condition is not as great as in the other test site 
conditions. Even though the BPF tone levels in the AWT are similar at Vq * 80 
knots and 180 knots, the higher harmonics of BPf are 2-5 dB lower at 180 knots 
relative to levels at 80 knots, the larger differences occurring at the higher 
tip speed Mach numbers. 

7.2 EFFECT OF BLADE NUMBER 


The spectra in Figures 40 through 43 also show that the higher harmonics 
of BPF decrease with increasing blade number, even though the effect of blade 
number is very small at BPF. The effect of blade number is shown explicitly 
in Figures 44 and 45. Figure 44 shows the sensitivity of BPF tone levels to 
Ng for the three test conditions, for the full-scale fan with low-flow 
operating line. Figure 45 shows the corresponding BPF harmonic spectra. 
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Figure 37. Component Source Contributions to BPF Tone for 
Full-Scale Fans in NASA-Ames 40 x 80 Wind Tunnel 
at 80 Knots. 
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Figure 39. Component Source Contributions to BPF Tone PWL for 
Full-Scale Fans in Flignt Environment at 180 Knots 
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Figure 45. Effect of Blade Number on BPF Harmonic Spectrum for 
Full-Scale Fan on Low-Flow Operating Line. 












Increaiing chc nunbcr of bUdct for a given solidity decreases the blade 
spacing and hence increases the tangential length scale-’to-blade spacing 
parancter lt;/st This results in a noise reduction^ as shown in Figure 11. 
However, increasing the nuisber of blades Hg also increases rotor blade 
aspect ratio h/c, resulting in a higher unsteady lift amplitude. The dipole 
noise source will therefore increase, off-setting the decrease due to increas- 
ing it/** dipole source only dominates at BPF (Figure 31), so the 

net result is little or no change in noise with Ng at BPF but decreasing 
noise with increasing Ng at second harmonic (n"2) and higher frequencies. 

7.3 EFFECT OF FAN SIZE 

To show explicitly the effect of fan size, some of the predictions were 
normalized with respect to fan rotor inlet area by subtracting 10 logj^O (A/Aj-gf) 
from predicted noise levels, %diere is a reference area, taken arbitrar- 

ily to be 1 m^. The BPF tone PWL results normalized in this fashion arc 
shown in Figure 46. The corresponding spectra are shown in Figure 47. The 
major observation to be made is that size does make a difference, i.c., rotor/ 
turbulence noise does not scale with fan area unless the turbulence scales are 
changed In proportion to the fan diameter change. 

The results shown in Figures 46 and 47 show that the smaller fans yield 
less noise than the area reduction effect can account for, by 5-10 dB. This 
again is related to the difference in tt/s which occurs because the turbulence 
scales (i^) remain the same while the geometric scales (s) decrease with de* 
creasing fan size. 

7.4 EFFECT OF OPERATING LINE 

Comparisons were made of rotor/ turbulence noise of full-scale fans with 
18 blades for two operating lines: (1) the low-flow (LF) line of Figure 28 

for typical high tip speed fans, and (2) the high-flow (HF) operating line of 

Figure 28 for typical subsonic tip speed fans such as the QCSEE fan. The BPF 

tone PWL comparisons are shown in Figure 48, while the spectral comparisons 

are shown in Figurv 49. It is seen from these results that, at a given tip 

speed Mach number, the HF fan is noisier than the LF fan. 
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Figure 46. Effect of Fan Size on BPF Tone PWL for 38-Blade Fans 
on Low-Flow Operating Line. 
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Figure 48. Effect of Operating Line on BPF Tone PWL for 
Full-Scale 18-Blade Fans. 
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The higher inlet relative Mach number for the high-flow fan is respon- 
sible for some of the noise differences shown in Figure 48, Assuming the 
noise to be proportional to a 3 dB higher noise level would be ex- 
pected for the HF fan based on Mach number level differences alone. The 
aspect ratio for the HF fan is higher because of the lower tip solidity, see 
Table 6, The unsteady lift is therefore higher for the HF fan. The aspect 
ratio effect alone is estimated to be 3.0 - 3.7 dB for this case. These 
effects account for most all of the differences shown in Figure 48 for the 
Outdoor Test Stand. 

For the Ames Wind Tunnel and Flight conditions, however, the higher axial 
Mach number for the HF fan gives a larger contraction ratio at a given tip 
speed. This has the effect of reducing the tangential length scale (t(./s), which 
increases the noise. However, it is also in the direction of reducing turbu- 
lence intensity at the fan face, which decreases the noise. The latter effect 
apparently overshadows the decreased-scale effect, since the HF fan and LF fan 
noise levels are nearly the same at low tip speeds. The rapid divergence near 
Mf; * 0.9 between the two operating- line curves is due to the fact that for 
the HF fan the rotor inlet relative Mach number is approaching unity at these 
high speeds, causing the quadrupole source contributions to become dispropor- 
tionately high as they approach their singularity point. 

The HF fan will produce more propulsive thrust at a given tip speed than 
the LF fan, because the flow is higher and the pressure ratio is the same. It 
is therefore of interest to compare the HF fan with the LF fan at the same 
thrust. For a typical LF fan approach power tip speed Mach number of M(- * 0.95, 
the static gross thrust (based on fan exhaust momentum) is 85,975 N (19,329 Ibf) 
The corresponding net thrust at 180 knots flight speed (gross thrust minus ram 
drag) is 51,735 N (11,631 Ibf). The equivalent thrust tip speeds for the HF 
fan occur at M^ * 0.78 and M^ = 0.815 for the static and 180 knots cases, 
respectively. These equivalent thrust points are Indicated on Figure 48 by 
closed symbols. It can be seen that the HF fan and LF fan rotor/turbulence 
noise levels are approximately the same at the same thrust. 
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8.0 NASA-AMES 40x80 FOOT WIND TUNNEL EVALUATION 


Based on the results of the parametric study discussed in the preceding 
section, an evaluation of the NASA-Ames 40x80 ft Wind Tunnel was carried out 
to assess the range of fan geometries, operating speeds, and tunnel speeds 
over which the Ames tunnel adequately simulates flight conditions as far as 
rotor/ turbulence noise is concerned. The objective was to define the oper- 
ational/ geometric boundaries for adequate flight simulation. 

Before carrying out the evaluation, an assessment of the wind-tunnel 
rotor/turbulence noise prediction level variability was made to establish 
an uncertainty band in predicted levels. As discussed in the preceding sec- 
tions, there is an uncertainty associated with taking turbulence properties 
measured 8-10 fan diameters upstream of the engine and projecting these to 
fan rotor inlet values using the sudden-contraction theory of Reference (5). 
One extreme is to assume that the sudden-contraction theory is correct, which 
qualitatively gives 

(ua/^a^fan ^'*a/'^o^up * C ^ 

and 

(ut/Ua)fan (ut/Vo>up * C-1/2 

Thus both axial and transverse intensities are lower after the contraction. 
This assumption gives the lower bound in noise predictions. The other extreme 
is to assume that the intensities do not change at all across the contraction 
(e.g., as in Reference 24), so that 

(^a/^a^fan “ ^“a/^o^up ^'^t/^a^fan * ^^t/'^o^up 

This gives an upper bound on the noise predictions. Now across a con- 
traction, U3 will decrease by some amount and u^. will increase. In any case, 
Ut will be substantially larger than u^, and hence the noise level will be 
primarily governed by the variations in U|- and variations in u^ will have 
little or no effect. An estimate of the difference between the upper and 
lower bound limits can therefore be made by noting (Reference 2) that the 
noise varies as the square of the turbulence intensity. Defining APWL as the 
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difference between upper and lower bound predictions, where subscripts A and 
B refer to precontraction and postcontraction values, we can assume that 

4PWL = - 20 logjo - 20 logic -m 

The numerator in Equation 26 represents the lower bound value of transverse 
intensity at the fan face given by sudden contraction theory, while the 
denominator represents the upper bound value corresponding to no change in 
intensity. Equation 26 was used to estimate the difference in calculated 
noise levels between the above two assumptions concerning contraction effects. 
Also, calculations were made over a range of tunnel speeds using the computer 
program directly; first with the sudden-contraction theory, and second with 
(ut/Ua)B * (^t/^o^A* difference between the two calculated levels was 

compared with the estimated difference given by Equation 26 above, and the 
agreement was within ±0.5 dB. 

Equation 26 is shown plotted in Figure 50, along with range of contrac- 
tion ratios encountered in this study for typical fan designs. It can be seen 
that the uncertainty band due to contraction effects is about 2-6 dB at Vg » 

80 knots over a range of tip speeds 0.5 £ Mj- ± 1.0. At 180 knots, the uncer- 
tainty band is only 0-2 dB for the tip speed range 0.75 ^ Mj- 1.0. For 
Mj; < 0.75 at Vg “ 180 knots, the contraction ratio is less than 1.0, i.e., 

Vo > Ug, and the applicability of contraction theory to expanding flows is 
questionable . 

It is emphasized that the uncertainty band given in Figure 50 is not due 
to variability and/or randomness in the tunnel turbulence properties, but is 
only due to the uncertainty in the techniques being used to extrapolate that 
turbulence information to the fan face values. In fact, the same uncertainty 
applies to the flight case, because upstream atmospheric turbulence charac- 
teristics are being extrapolated to fan face values with the same techniques. 
Also the maximum (6 dB) uncertainty band is less than the tone fluctuation 
amplitudes usually observed for rotor/turbulence noise (References 21 and 22). 
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The outdoor test stand predictions also have an uncertainty band due to 
the variability in turbulence properties which are measured at the fan face. 
For example, measured axial length scales varied from 86 to Z59 ft (26 to 
79 m) at the GE Peebles, Ohio test facility for one engine test. A parametric 
calculation of the influence of An was made (An/A^ held constant) for the 
full-’Bcale low-flow fan, and these results are shown in Figure 51. From 
these results, it can be concluded that an uncertainty of tl dB to -4 dB can 
exist in predicted Outdoor Test Stand levels due to uncertainty in knowing 
axial length scale, recalling that £ nominal value of A^ *■ 100 ft (30.5 m) 
was used for all OTS calcuUtions discussed in the previous section. Also, 
the measured transverse intensity u^/Ua varied by a factor of two, implying 
an additional t3 dB uncertainty due to variability in U{;/Ua. 

To illustrate the above effects, measured values of BPF tone forward- 
arc PVn. on a CF6-50 engine were compared with predicted rotor/turbulence 
noise tone levels using the above uncertainty estimates to define a band 
of expected levels. This comparison is shotm in Figure 52. Even though 
the engine had an acoustic suppression liner in the inlet duct whose effect 
on rotor /turbulence interaction tones is quantitatively unknown, the data 
is seen to fall within the prediction band. 

Finally, the uncertainty in predicted flight rotor/turbulence noise 
levels must be considered. Houbolt, Reference 25, quotes many data sources, 
and suggests that a rather alarming variability in atmospheric turbulence con- 
ditions can prevail, depending upon the proximity of weather fronts, wind 
shear, squall lines, thunderstorms, temperature gradients, 'terrain, etc. Air- 
craft wake turbulence in the vicinity of airports will introduce variations 
dependent upon traffic patterns, density and duration. Houbolt estimates 
the variation in atmospheric turbulence length scales to be 300 to 700 ft (91 
to 213 m). A nominal value of 500 ft (152 m) was used in, the present study. 

A parametric calculation of the effect of length scale (precontraction) is 
shown in Figure 53 for a full-scale, high-flow (18 blades) fan. Again it is 
seen that the uncertainty in predicted noise levels for the length scale range 
of 300 to 700 ft is about ±2 dB. Houbolt also states that atmospheric turbu- 
lence rms velocities of 3.0 - 3.5 fps are average values, with possible peak 
"gust" values of 20 times as much. Such extremes are rare and are confined 


97 








watts 
















to 80 vure weather conditions, but nuvertheluss highlight the variability that 
can bo expected in rotor/turbu^ence noise in flight. A 4:1 variation In pre~ 
contraction uj, is not unreasonable in flight, which would then yield a ±6 dB 
variation in noise. Combined with the i2 dB uncertainty due to variations in 
length scale, a total uncertainty band on predicted flight rotor/turbulence 
noise levels of ±8 dB is the best that can be expected, based on atmospheric 
turbulence variability from the "expected average" values. 

From the results shown in Figures 44 through 49, it is seen that rotor/ 
turbulence noise levels in the Ames 40x80 Wind Tunnel are predicted to be 
15-20 dB lower, on the average, than the Outdoor Test Stand levels, but still 
about 10-15 dB higher than corresponding Flight levels. The question that 
arises is, are the Wind tunnel levels low enough? To answer this, first 
recall that the noise level increases with decreasing blade number, and that 
full-scale fans produce more "noise-per-unit-area" than do scale model fans. 

If we take into account all of the above uncertainty band limits and consider 
as the worst case the full-scale 18-blade fan, the comparison of predicted 
BPF tone levels for the various test conditions would look more like Figure 54, 
where a band or range of values is shown rather than a single average line. 

The point to be learned from this figure is that when variability is taken 
into account, the wind tunnel levels can bo as low as flight levels, even 
though the average line predictions show a 10-15 dB difference. Also, the 
uncertainty band associated with variability in turbulence properties for 
the Outdoor Test Stand and Flight conditions is considerably larger than the 
uncertainty for the Ames 40x80 Wind Tunnel due to extrapolating upstream turb- 
ulence characteristics to the fan face conditions. 

Another way of assessing the adequacy of the NASA-Ames 40x80 Tunnel for 
simulating flight turbulence levels is to see how much tone reduction is 
obtained relative to the fan characteristic broadband level. An empirical 
method for correlating fan broadband noise was published in Reference 28, 
and this method was used to correlate the broadband levels of several scale 
model fan stages, using data with a turbulence control structure in place 
(Reference 6). The correlations, given in Reference 29, relate the peak 
broadband noise level to rotor inlet tip relative Mach number Mj. and rotor 
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tip incidence angle (relative air angle minus blade leading edge cainberline 
angle). This correlation was used to estimate the broadband level (1/3 octave) 
at the blade-passing frequency, and these results arc also shown in Figure 54, 
It can be seen that the NASA-Ames 40x80 Wind Tunnel Rotor/Turbulencc tones 
arc well below the predicted broadband levels at all but the highest tip 
speeds. 

The spectral distribution of the rotor/turbulence tones is shown in 
Figure 55, compared with predicted l/3-octave fan broadband noise spectra. 

It is apparent from this comparison that the higher harmonics of BPF rotor/ 
turbulence tones are well below the fan broadband noise level for the Ames 
40x80 Wind Tunnel case. Additional higher harmonic tone sources (i.e., 
rotor-stator interaction noise) will no doubt clominate the spectrum. 

Finally, a check was made on the ability of the NASA-Ames 40x80 wind 
tunnel to simulate an aircraft accelerating taxi or "ground roll", where fan 
speed is held constant and wind speed or aircraft speed is increased from 
zero to some typical approach speed. Predictions were made of tone PWL 
for the JT15D engine (see DATA/THEORY COMPARISONS section) at a fan tip speed 
of Mt^"0.93, over a range of tunnel velocities from 10 to 200 knots, using the 
sudden contraction theory (Reference 5). Taking into account the possibility 
that the transverse turbulence intv nsity at the fan face does not follow the 
sudden contraction theory but instead remains constant (Reference 24), addi- 
tional calculations were made using Equation 26. 

A "ground roll" calculation was also made, by specifying an initial 
precontraction static turbulence characteristic. This characteristic was 
chosen to yield approximate Outdoor Test Stand fan face turbulence at a 
very large contraction ratio. The flight speed was then progressively 
increased, causing contraction ratio to decrease, all the while holding pre- 
contraction turbulence properties constant. 

The precontraction turbulence properties were selected based on the fol- 
lowing reasoning. In an outdoor test stand, the transverse intensity at the 
fan face is typically 4%, which, for Ug s 400 fps at approach power setting, 
gives a fan face turbulence rms velocity of ~16 fps, From the data given by 
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Houbolt (Reference 25), atmct-oheric turbulence levels are about 3 fps on a 
normal day. The implied transv.^rse velocity ratio is 16/3 or 5.33, and sudden 
contraction theory then suggests that, for this velocity ratio, the contrac- 
tion ratio is about C z 40. Assuming a typical measured axial length scale 
inside the fan duct of -100 ft, the corresponding transverse scale is -0.4 
ft, and the scale ratio 1^1 - 250, about what has been measured in Refer- 

ence 6. The implied atmospheric length scale is g, z 2.5 ft, not unrealistic 
for ground level turbulence. The sudden contraction theory then gives an axial 
turbulence intensity at the fan face of -0.00053, much lower than the usually 
measured 0.005 to 0.010. However, it is the transverse intensity at the fan 
face which determines the noise level, so the underestimation of the axial 
intensity is of little consequence. It is also possible that these lower 
values of axial intensity cannot be properly measured due to instrumentation 
noise floors, contamination by acoustic velocity fluctuations, etc. 

Figure 56 shows the results of the calculations described above. The 
"ground roll" curve shows a rapid decrease in tone level from the Static 
level, reaching about 17 dB reduction at -100 knots taxi speed. A similar 
effect was measured by Lowric (Reference 30), his data showing approximately 
15 dB reduction in BPF tone level at 80 knots, in an accelerating taxi test of 
a VC10-RB211 flying test bed. NASA Ames 40x80 Wind Tunnel test data on the 
JT15D engine (Contract NAS2-8675 with GE, unpublished data) is also shown in 
Figure 56. Conside 7 fing the possible variability in atmospheric turbulence 
properties and the reasonable agreement between "ground-roll" predictions 
and both predicted and measured wind tunnel levels shown in Figure 56, it 
can be concluded that the 40x80 wind tunnel is an adequate simulation of 
rotor/turbulence noise forward-speed effects. 

In summary, based on all the calculations and parametric studies dis- 
cussed herein, it is concluded that the NASA-Ames 40x80 Wind Tunnel provides 
adequate suppression (or reduction) of rotor/turbulence noise such that proper 
simulation of flight fan source noise characteristics is obtained, for sub- 
sonic tip speeds. Near Mj. approaching unity, wind tunnel rotor/turbulence 
tones may be as high as or higher than the broadband level, but cut-on of the 
rotor-alone noise field at Mj. = 1 will probably mask this effect. 
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Figure 56. Comparison of Predicted "Ground Roll" Rotor/ 
Turbulence BPP Tone PWL with Predicted and 
Measured Wind Tunnel Turbulence Levels. 
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9.0 CONCLUSIONS 


The following conclusions were drawn relative to the results obtained 
from this study: 

1. The NASA-Ames 40x80 Wind Tunnel adequately simulates forward velocity 
effects on fan noise inlet arc tones produced by rotor/turbulence 
interaction, for fans from 0.5 to 2.0 m in diameter, for all subsonic 
tip speeds, and for blade numbers from 18 to 38. The simulation 
should be adequate for higher and lower values of diameter and blade 
number, but the present study was confined to the above ranges. The 
wind tunnel rotor-turbulence interaction tones are sufficiently low 
in level such that internally generated fan noise sources dominate 
the fan noise observed spectrum, and the rotor-turbulence noise has 
no contribution to the observed spectrum. 

2. The adequacy of utilizing sudden contraction theory to predict fan 
inlet turbulence spectra from measured turbulence properties at some 
distance (several diameters) upstream of the fan inlet is question- 
able. 

3. Rotor/turbulence interaction noise produced by ingestion of large- 
scale turbulence in a contracting inflow is primarily a narrow-band 
tone source, with no significant contribution to the fan broadband 
level. 

4. For large scale inlet turbulence (i>D|;), the associated noise levels 
are primarily determined by the ratio of transverse turbulence scale 
to blade spacing at the fan face, and are not materially affected by 
axial scale if transverse scale is maintained constant. Maximum 
noise occurs when t^/s<l. 

5. Rotor/turbulence noise in Che NASA-Ames 40x80 Wind Tunnel is 15-20 
dB below that produced in an outdoor static test stand for wind tun- 
nel velocities greater Chan 40 kCs. The reduction is due Co the 
much smaller turbulence scales and intensities in Che wind tunnel 
environment. 

6. The Rotor/turbulence noise In Che NASA-Ames 40x80 Wind Tunnel is 10^' 

15 dB higher Chan Chat produced in flight at the same flight speed, 
but still several dB lower Chan Che other sources of fan noise, e.g. 
broadband, rotor-alone multiple-pure-tones, rotor-stator interaction, 
etc. Even though Che flight atmospheric turbulence intensities are 
higher, Che turbulence scales are so large at altitude that 

Che short wavelength components of Che turbulence spectra have very 
low turbulence energy and so yield correspondingly low noise levels. 

7. Under contracting inflow conditions (Ua>V(,), the transverse tur- 
bulence intensity controls the level of generated rotor/turbulence 
noise, and Che axial intensity has a relatively small influence. 
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A rather wide uncertainty band in expected rotor/turbulence noise 
levels exists, :f5~8 dB, for outdoor static and flight conditions, 
due to the random variability in turbulence conditions. The wind 
tunnel turbulence characteristics are, on the other hand, much more 
consisv';ent. 
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APPENDIX A 


LONG AXIAL LENGTH SCALE APPROXIMATION 


As ditcuBsed in Reference 2, the rotor/ turbulence noise spectrum theo- 
retical model employed in the present study consists of a summation of compo- 
nent spectra centered at blade-passing frequency and its harmonics, n “ 
1,2,3...,., etc. The component spectra have shapes which ore essentially the 
turbulence spectrum shape as axial wave number k| is varied over the range 
- • j< 1*1 £ ♦ •» with integrals over lc 2 and k 3 performed at each value 

of k^. This component summation is illustrated qualitatively in Figure 57. 

# 

For very large length scales i^/s » 1, the component spectra 1-1, 2-2, 

3-3, etc., shown in Figure 57 become very "peaky" and narrow, and the contri- 
butions of component spectra adjacent to a spectrum centered at f"nf|, are negli- 
gible in the vicinity of f"nf|,. For example, the 1-1 and 3-3 component 
spectra in Figure 57 would contribute a negligible amount to the total spectrum 
at and around f*2f|,, and the 2-2 component is all that is required to evaluate 
the total spectrum level around f^lfi,. 

When the spectrum around f*nf|} is very narrow, as occurs when i^/b » 1, 
the total power level contained in the peak can be evaluated by integrating 
over f and hence kj. If the turbulence spectrum drops off very rapidly with 
k]^, compared to variations in either unsteady lift amplitude or quadrupole 
source amplitude, these factors can be approximated by their values at kj^ « 0 
and the integral of only the turbulence spectrum over kj^ need be taken. 

Further, the spectrum levels need only be evaluated at harmonics of blade- 
passing frequency f“nfj,, where n * I, 2, 3, ....etc., and only the component 
spectrum centered at f^nf],, i.e., the n-n component, need be calculated. Thus 
no summation of component spectra is required. 

The integral of the turbulence spectra over kj^ can be derived in closed 
form from equations (16). The resulting expressions are as follows: . 
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♦Il<V - I ^1 


♦ j2 ^•‘2^ “ y* ^12 ^*‘r **2^ d kj - 0 


♦ 22 <>^ 2 ^ 


y* ^22 


k2> d kj 
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t t 
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ita: 


where 


1 ♦ 4 


(27) 


(28) 


(29) 


(30) 


Thus, in the long length-scale approxination » 1, the total sound 

power in the vicinity of f <■ nfb is calculated by evaluating the unsteady 
lift dipole and quadrupole source strengths at k|^«0 and multiplying by a 
suitably weighted combination of the turbulence spectra given by equations 
(27-30) above. 

It can be seen from equations (27-30) that the turbulence spectra for 
f>g/B l,i.e., integrated over k^, are independent of This is con- 

sistent with the results shown in Figure 10 which indicate that the noise 
level approaches a constant value independent of A^/s for A^/g » 1. This 
occurs because the peak value at f « nf|, (Figures 4 and 7) varies linearly 
with A^/s for Ag/s » 1, &i;ereas the effective bandwidth decreases inversely 
with Ag/s. 

VAien the turbulence undergoes a contraction, the long length-scale 
approximation can still be employed. Referring to equations (19-21), the 
postcontraction spectra and 4>22(JS) only need to be considered since 

for A/s >> 1 the Integral over kj_ (or Kj^) of * 12 (jl^ vanishes. Now the 
Integral over of and ^22 Involves the expression 


115 


f <> <1 - J f <1 kl 

and a •imilar one for # 22 * examination of equations 

(19-21) that only the leading terms will contribute significantly since they 
involve the integrals /^n dki and /♦22 ‘^•^1* whereas the remaining terms 
involve / which should be small if fij drops off rapidly enough 

with ki. The postcontraction spectra corresponding to equations (27-29) can 


therefore be approximated by the following: 

00 

♦11 («2) " JJ *11 dK3dKi 

^ ♦u 


or 


<*^ 2 > ^ 

Similarly, 


♦11 (k 2 ) 


hz ^ ^ (^ 2 ) 


(31) 


(32) 


where k 2 ■ < 2 /^* as before. Equations (31,32) are used In place of (27,29) 
to account for contraction effects for C ■ Ug/Vo ^ !• 
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APPENDIX B 


SYMBOLS AND ABBREVIATIONS 


A - Fan rotor inlet annulua area, alao turbulence apectrum variable 

AWTl - AaK!a wind tunnel, Vg ■ 80 knots 

AWT2 - Aiaea wind tunnel, V© ■ 180 knots 

a - Conatant 

BPF - Blade-paaaing frequency 

C - Contraction ratio 

C| " Rotor blade lift coefficient 

CR - Contraction ratio 

c ~ Rotor blade chord 

DV - Diacharge valve (throttle) aetting 

D(; - Fan rotor tip diameter 

F(k) • Turbulence apectrum function 

FLT - Flight, Vo - ISO knots 

f • Frequency, alao isotropic turbulence correlation function 

fg > Blade-passing frequency 

C(k) - Turbulence spectrum function 

g - Isotropic turbulence correlation function 

HF - High-flow operating line of Mg versus 

HS - Half-scale 

h - Blade height 

k - Precontraction wavenumber 

-mm- 

k£ - i*”^ component of k 

L - Turbulence integral length scale 

LF ■* Low flow operating line of Mg versus 

i - Isotropic turbulence length scale 

Ag - Axisymmetric turbulence axial length scale 

At - Axisymmetric turbulence transverse length scale 

Mg - Fan inlet axial Mach number 

Mt - Fan inlet tip speed Mach number ' 

Mf - Fan inlet relative Mach number 
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Nb 

Np 

n 

OTS 

P 

Po 

Pi 

Pt 

PWL 

PWLo 

PWLi 

PWLt 

APWL 

Q1»Q2»Q3 

QS 


rt 

SAC 

s 

TCS 

U 

Ua 

Ut 

ua 

“t 

Vo 

*i 

i’ij 

6ij 


rotor blade number 

fan rotative speed, rpn 

blade-passing frequency harmonic number 

outdoor test stand 

acoustic power, watts 

peak acoustic power, watts 

acoustic power at band edge, watts 

total acoustic power in band, watts 

acoustic power level, dB re: 10**^^ watts 

decibel equivalent of Pq 

decibel equivalent of 

decibel equivalent of 

PWLo - PWLl 

axisyinmetric turbulence correlation functions 
axisymmetric turbulence correlation tensor 
quarter-scale 

turbulence velocity correlation function 

spanwise radial coordinate 

rotor inlet hub radius 

rotor inlet tip radius 

rotor inlet rms pitchline radius 

Schenectady anechoic chamber 

blade-to-blade spacing 2irr/Ng 

turbulence control structure 

local mean flow velocity 

fan rotor inlet axial velocity 

fan rotor inlet tip speed 

axial component of rms turbulence velocity 

transverse component of rms turbulence velocity 

wind t'znnel velocity or flight speed 

£th component of separation vector x 

transverse separation 

turbulence spectrum coefficient 

kronecker delta 
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*ijk 

e 

K 

\ 

U 

T 


♦ij 

♦ij 



> alternatiing unit tensor 

- integration variable 

- poatcontractlon wavenumber 

- i^h component of ic 

- length scale ratio t^/l^ 

“ turbulence velocity ratio U|-/ua 

- cross-correlation time delay 

- precontraction three-dimensional turbulence spectrum 

- precontraction two-dimensional turbulence spectrum 

- postcontraction three-dimensional turbulence spectrum 

- postcontraction two-dimensional turbulence spectrum 


Subscripts 


1 - axial direction 

2 - tangential direction 

3 - radial direction 
0 - free stream 

a - axial 

t - transverse, tangential^, tip 
b - blade 
p - pitchline 
h - hub 
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